Mode coupling induced dissipative and thermal effects at long times after a quantum 

quench 
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An interaction quench in a Luttinger liquid can drive it into an athermal steady state. We 
analyze the effects on such an out of equilibrium state of a mode coupling term due to a periodic 
potential. Employing a perturbative renormalization group approach we show that even when the 
periodic potential is an irrelevant perturbation in equilibrium, it has important consequences on 
the athermal steady state as it generates a temperature as well as a dissipation and hence a finite 
life-time for the bosonic modes. 
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The high degree of tunability and control associated 
with cold-atomic gases make them an exciting testbed 
for studying a host of phenomena related to interacting 
quantum particles Among the topics of great cur- 
rent interest are the non-equilibrium physics of quantum 
quenches for which these systems are particularly well 
adapted (see [3(] and references therein) and other classes 
of steady state non-equilibrium phenomena such as sys- 
tems subjected to a time dependent noise 

Crucial questions in all these out of equilibrium phe- 
nomena is the relaxation mechanism by which the system 
reaches a steady state and the properties of the steady 
state, in particular whether the latter is athermal or ther- 
mal and hence described by a Gibbs distribution. In the 
case of a quench where non-trivial time-evolution is trig- 
gered after a sudden change in system parameters, many 
models mainly in one-dimension [ij and involving sim- 
ple effective theories @, @] , are found to reach an ather- 
mal steady state characterized by a generalized Gibbs 
ensemble (GGE). However the generality of the GGE re- 
mains under debate since certain observables do not obey 
it More generally, the relaxation mechanism 

and the nature of the steady state in more complicated 
field-theories, as well as the role of a finite system size in 
numerical studies, are still largely unknown 
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It is thus important to have theoretical models for 
which such non-equilibrium questions can be reliably 
studied. One good candidate for such an analysis is a one 
dimensional system of interacting bosons, leading to the 
so called Luttinger liquid physics The excitations 

of such a system can be represented by density modes 
which are essentially independent. On such a system, 
quenches corresponding to a change of the interaction 
reveal a steady state which still has independent modes, 
but these are now characterized by a non-equilibrium dis- 
tribution that does not relax to a thermal state • 

In this paper we examine the effect of a mode coupling 
term on the above system, where the mode coupling is 
due to an externally imposed optical lattice. We address 
explicitly the question of thermalization and asymptotic 



relaxation of such a system. We use the Keldysh tech- 
nique L^ and a controlled renormalization analysis and 



show that even in cases for which the lattice potential 
would be irrelevant in equilibrium, it leads for the out of 
equilibrium situation to the appearance of finite dissipa- 
tion, as well as a finite temperature for the low energy 
modes. We thus explicitly obtain a mechanism, which 
we argue should be generic, in which thermalization and 
dissipation arises due to the transfer of energy from the 
long wavelength modes to the short wavelength modes, 
the latter thus acting as a bath. 

We consider interacting one dimensional bosons in the 
continuum. The low energy properties of such a system 
can be efficiently represented by a Luttinger liquid [l3| 
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where <j> is related to the long wavelength part of the 
density by p{x)—— V4>(x)/tt, while LI is the canonically 
conjugate variable to <f>. The eigenmodes of the Hamil- 
tonian are the sound waves of density with a dispersion 
oj=uoq. The information about the interactions is con- 
tained in the two Luttinger parameters: uo the velocity 
of density oscillations, and Kq a dimensionless parameter 
controlling the decay of correlation functions. 

The bosons are driven out of equilibrium via a sudden 
interaction quench which for the effective Luttinger liquid 
theory, simply implies a sudden change of the Luttinger 
parameter from Kq — > K, and the velocity from uq — > u. 
To satisfy Galilean invariance we choose u=vpjK and 
Uo=vp/ Kq. The time evolution of the initial state is 
therefore due to Hf=Hi(Ko — > K, uq — s- u). 

We first give here the full solution for a quench Kq — > 
K with arbitrary interactions. In the Keldysh formal- 
ism it is convenient to define classical {4> c i — {4>- + 
4> + )/\/2) and quantum ((^> g =(</>_ — ^ + )/v2) fields where 
4>-(+) are the time (anti-time) ordered fields on the 
Keldysh contour. In terms of these fields, the action 
that describes the steady state behavior at long times 
(t + t' — > oo ) after the quench when transients related to 
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FIG. 1. C K {q = 0.5, cj) and 2Im[C H (g = 0.5, u)] in equilib- 
rium K eq = if,ie g = 2.0 and after a quench K eq — 2, -ffneg = 3 
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Eq. ([2]) implies that the retarded propagator — i{ 
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— . r is identical to that in the 

(cj+zo) — u A q A 

ground state of Hf, while the Keldysh propagator G K = 
—i(<j) c i<j)* cl ) which is sensitive to the occupation of the 
bosonic modes is, 
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G (q, oj) = — ^ 1 + I sign( W ) [G R - G A ) (3) 
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Thus the fluctuation-dissipation theorem (FDT) defined 
by G K = (G R - G A ) coth(cj/2T), where T is the tempera- 
ture of the bosons, is violated. When K=Kq, FDT is re- 

T— 

covered as coth(w/2T) > sign(cj). Note that although 

the system is now out of equilibrium, each q mode is still 
infinitely long lived since <5=0 + . 

Sb can be used to evaluate the basic two- 
point correlation functions corresponding to the den- 
sity fluctuation Cj^ and response defined 

^((0_(r,t)-0 + (O,O)) 2 ) 



as Cfo(r,t) = 
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e -i-<(0_(r,t)-0 + (O,O)) 2 ) 

K (K 2 /K 2 + 1), we find, 

Xl<(0_(r,t)-0 + (O,O)) 2 ) = 
^lnf ° 2 + <".; + " )2 V^lni 
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where a is a short distance cut-off, and 7 is an arbitrary 
coefficient which will later be related to the periodicity 
of an externally imposed lattice potential. In equilib- 
rium K eq =K neq , and one recovers the usual power-law 
decay of Luttinger liquids. However out of equilibrium 
one finds power-law behavior with new decay exponents 
K neq . For the case of Kq=1 this power- law decay was ob- 
tained in [1| , and can also be recovered using a GGE that 
accounts for the conservation of the occupation number 
of appropriate bosonic modes. Since K neq > K eq , the 
propagators always decay faster than in equilibrium. 

The role played by the oscillating terms in Eq. (HJ 
which differentiates between response and correlation 
functions has not been explored before, and will play an 
important role in the RG. Its importance can already be 
seen at this level by studying the FDT ratio defined by 
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this ratio reduces to 



the equilibrium T=0 result of sign(w), out of equilibrium 
it can be used to formally define a oj,\q\ dependent ef- 
fective "temperature". In the limit uj —> 0,q —} 0, the 
effective temperature T defined as 
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and therefore assumed to be T > uj is, 
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where the energy-scales are expressed in units of u/a, 
and length scales in units of a. The behavior of 
C K (q=0.5,uj) and 21m Cf^q = 0.5, oj) are plotted in 
Fig. [T]for the equilibrium case when K eq =K neq =2 and 
the non-equilibrium case of K eq =2, K neq =3. Note that 
this "temperature" is dependent on the correlation func- 
tion we use, contrary to the case of equilibrium for which 
each ratio between response and fluctuation defines the 
same temperature. Fig.[Tjshows that besides the appear- 
ance of an effective temperature, a striking effect is the 
appearance of a dissipation characterized by a non-zero 
slope of Im.[C fl ] oc —ir\uj. As we shall show below, the 
temperature and the dissipation already apparent at this 
stage will reappear in the RG analysis. 

We now study how this non-equilibrium state is mod- 
ified by a coupling between modes. Although in prin- 
ciple any form of nonlinear coupling, such as for exam- 
ple a (j) 4 term can be used, we focus here on the case 
of a cos(70) perturbation. There are two reasons for 
such a choice: i) If the phase 4> represents real inter- 
acting bosons, the Hamiltonian cannot contain pertur- 
bations coupled directly to <j> but only to derivative or 
periodic functions of 4>] ii) such a periodic term arises 
naturally when a periodic potential is added on the sys- 
tem. It is the source of the Mott transition in one dimen- 
sion [l3j |. and thus very natural to study in that context. 
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As for the case of equilibrium we study this term by a 
renormalization group (RG) procedure, since the rest of 
the Keldysh action S is quadratic. The Keldysh path- 
integral is Zx = jT> [(p c i, 4> q ] e i ( So+Se 9> where 



x10 



dx / dt [cos (70-) — cos(70 + )] (7) 



a _9^_ 

dsn — n 

a 2 

Writing such an action assumes that after the quench 
has long taken place, one switches on the cosine term 
infinitely slowly. Without the quench the system would 
thus relax to the ground state (at T=0) in the presence 
of the cosine term. To which state the system will tend if 
the initial state is not in equilibrium is the very question 
we address here. In order to perform an RG analysis, 
we split the modes between slow and fast components 

(t>cl,q( xt ) = <t>cl,q{ xt ) + 4>cl,q( xt ) ^eie 
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and A/A' = e d ( lnL \ and integrate out the fast modes. 
Care has to be taken to regularize the two-point functions 
appropriately, and no cutoff on time should be imposed. 
We choose here a standard momentum cutoff [l3[, the 
details of the computation will be given elsewhere. 

The first step of RG generates a correction to S sg as 
well as corrections to So which may be absorbed into a 
redefinition of K and u. However when K ^ Kq, addi- 
tional terms of the form 
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are generated. These corrections can be summarized in 
the following RG equations, 
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FIG. 2. Strength of the dissipation n for _Ko=3, 7=2 and 
g=0.05 and 3=0.1. 



lK/u=fZc dr C dt(r 2 Tt 2 )Im e -V((*-(*.'-)-*+(o,o))' , >j? 
and F=\ (e lA (*+'') + e lA ^) + 

{sin(A(i + r)) + sin(A(t - r))} 



Eq. (flOl) reflects the new scaling dimension of the op- 
erator e %1 ^ due to the change in the decay exponent from 
K eq to K neq . It still defines two regimes, one in which the 
cosine is irrelevant, and one for which the perturbative 
RG would lead to strong coupling for the cosine term. In 
equilibrium this would reflect the Berezinski-Kosterlitz- 
Thouless transition corresponding to the Mott transition 
(K = 2 and infinitesimal g). Eq. (ITTI) is the usual scal- 
ing of K, which is reduced by the presence of the cosine 
term. This equation is also slightly modified compared 
to the equilibrium one when K eq 7^ K neq . Eq. ([T2"j) is 
a renormalization of the velocity. It appears here be- 
cause we took a pure momentum cutoff which thus does 
not respect the Lorentz invariance. It would appear also 
in equilibrium with the same cutoff structure. These 
three equations would thus lead to two separate phases, 
one in which the cosine is irrelevant, and a strong cou- 
pling regime whose physics would be beyond the reach 
of the perturbative RG. In order to stay in the regime 
for which the RG is reliable even asymptotically we con- 
centrate here on the case K eq , K neq > 2 for which the 
cosine term is irrelevant according to Eq. flTU|) . Other 
regimes of the phase diagram will be discussed elsewhere. 
In this regime one could naively expect to recover the 
same physics as without the cosine (namely the ather- 
mal state corresponding to Eq. Q). However the two 
remaining equations (|13I14[) introduce qualitatively new 
physics and lead to quite a different state. 

Eq. (fTBf shows that contrary to the case of an equilib- 
rium quantum system, for which the friction coefficient 
remains always infinitesimal (77 — + ), even at finite tem- 
perature, here because of the combination of the cosine 
term and the initial out of equilibrium action, a finite 
friction is generated. If one starts from the equilibrium 
situation K=Kq then of course /t,»7=0 and one recov- 
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ers the conventional results. The finite friction causes a 
crossover of the mode dispersion at low energies from a 
pure quantum behavior, dominated by (d t (f>)' 



to a more classical one r]dt4> ~~ ^ inuxj), and the correlation 
functions will reflect this. Interestingly the physics of 
dissipation can also be recovered in a quench involving 
fermions [15j . Similar to the case studied here, an ini- 
tial quench on a system of noninteracting fermions can 
cause it to reach a non-equilibrium steady state charac- 
terized by a highly broadened distribution function. As a 
result switching on infinitesimal interactions (which can 
be treated within the random phase approximation) can 
cause efficient scattering and an enhanced particle-hole 
continuum which leads to a damping of collective modes, 
at least for attractive interactions between fermions. 



In addition to the generation of the friction, Eq. (|14j) 
shows that a constant term (in the limit cj — > 0) 
is added to the Keldysh part of the action. In this 
limit and in equilibrium, this term would simply be 

oc uorj coth "~°> 2Tr]. Thus the constant term to- 
gether with the appearance of a finite friction can be 
interpreted as a finite temperature, at least for small 
enough frequencies. One thus recovers at small frequen- 
cies the action (so called Martin-Siggia-Rose action) of a 
classical system with a finite friction and a thermal noise. 
Note that because we have shown that the full action 
renormalizes to a quadratic one, and that the Keldysh 
term tends to a constant, the temperature as defined 
above is indeed the one that will appear in all 4> corre- 
lations, at least asymptotically for low frequencies, con- 
trary to the case of ([2]) . Therefore the non- linear coupling 
of the modes leads to a thermalization of the system. 
The full frequency dependence of such a noise is however 
quite complicated leading to an interesting crossover, de- 
pending on the frequency scale, between the athermal 
distribution and the classical, finite temperature one at 
low frequencies. In particular the RG flow itself has been 
derived from the quantum athermal correlations. The 
corrections generated by the RG will thus change signifi- 
cantly at a scale for which uj 2 = , q{uj)uj. Since at this scale 
the system enters in a more classical regime with expo- 
nentially decaying correlation functions (see below), this 
regime will not change the fact that the cosine is irrele- 
vant, and will simply slightly modify the final values of 
the friction and temperature. 

Fig. [2] shows the solution for the renormalized 77 for 
two different g and Kq=3. The non-monotonic behav- 
ior arises because the larger is K the more rapidly g 
renormalizes to zero leading to a smaller renormalized 
77. While at the other end, when K=K , rj=0. These 
two behaviors have to go through a maxima. Quite nat- 
urally the friction is proportional to g 2 . This is however 
not the case for the temperature T e ff which for small 
g is found to reach the following value independent of g 



(where we have set F=l in the RG equations) 
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While as K neq K eq , rj,nT e ff 



0, T e ff is non-zero 
This is because the order of limits w — > 0, K eq — > K n 



eq 



do not commute. Further T*j^K* eq /K* q ~ T, and hence 
is consistent with the non-interacting estimate for the 
temperature ([6]). 

Let us finally compute the correlations at the thermal 
fixed point where the action is (dropping ui 2 terms in 
comparison to ur/) 
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The above implies that equal-time two point correlation 
functions decay exponentially in position, 
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while the dissipation affects unequal-time correla- 
tion functions G R (q,t)=~6(t){nK*u/ri*)e- u2 <? t ^\ 
G K (q,t) = _(2«ir)(ZM^) e -«V|t|/')\ Thus in an 

experiment involving a one dimensional Bose gas in a 
periodic potential [l|, a probe of the density-density 
response function which directly correspond to the 
correlators (e^e -17 ^), should reveal the dissipative and 
thermal effects we predict. 

In summary, by studying the particular example of a 
quenched Luttinger liquid in the presence of a lattice, we 
have found a mechanism, that we believe is generic by 
which a non-equilibrium system in the presence of mode 
coupling will both thermalize and acquire a finite fric- 
tion or lifetime for the modes. It is important to note 
that these effects are related to the presence of a contin- 
uum of excitations in the system by which local degrees 
of freedom can relax and exchange energy. By this argu- 
ment it is possible that thermalization might not occur 
in the Mott insulator phase, and the fact that T e ff van- 
ishes near the critical point, might be a prelude to this 
physics. An investigation of this issue and also how the 
results depend upon the rapidity with which the cosine 
potential is switched on, are important open questions 
left for future studies. 
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